A numerical algorithm for MHD of free surface flows at low magnetic Reynolds numbers

We have developed a numerical algorithm and computational software for the study of magnetohydrodynamics (MHD) of free surface flows at low magnetic Reynolds numbers. The governing system of equations is a coupled hyperbolic/elliptic system in moving and geometrically complex domains. The numerical algorithm employs the method of front tracking for material interfaces, high resolution hyperbolic solvers, and the embedded boundary method for the elliptic problem in complex domains. The numerical algorithm has been implemented as an MHD extension of FronTier, a hydrodynamic code with free interface support. The code is applicable for numerical simulations of free surface conductive liquids or flows of weakly ionized plasmas. Numerical simulations of the Muon Collider/Neutrino Factory target have been discussed.